Quantum Phases of Ultra Cold Bosons in Incommensurate ID Optical Lattices 
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We study theoretically a BEC loaded into an optical lattice in the tight-binding regime, with a 
second, weak incommensurate lattice acting as a perturbation. We find, using direct diagonalization 
of small systems and a large scale, number conserving Gutzwiller-based mean field approach, the 
phase diagram of the model. We show that the superfluid, Mott-insulator and Bose-glass phases 
familiar from disordered systems appear here as well, with a localization transition occurring al- 
ready in a ID geometry, contrary to the well known Anderson localization. We complement these 
results with GPE simulations, and stress the fundamental difference between commensurate and 
incommensurate lattices. These results, together with calculations we have performed for realistic 
scenarios, suggest that this setup is well suited for available experimental realizations. 
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The possibility of observing quantum phase transitions 
by loading a Bose-Einstein Condensate (BEC) into an 
optical lattice has focused attention to disorder in- 
duced effects, based on the non-interacting Anderson 
phase transition Q • Although Anderson localization was 
observed (e.g. in microwave systems, see Q), the local- 
ization transition was not seen. Following the compre- 
hensive treatment of the interacting case in Q , research 
specific to BECs in optical lattices in the presence of dis- 
order has been done (e.g. 0, Q)- It was shown that 
the phase diagram of the system includes a compressible 
superfluid phase with a gapless excitation spectrum, a 
compressible (gapless) glass phase with no superfluid re- 
sponse, and an incompressible (gapped) insulator phase. 

Experimentally, attempts to create this form of ran- 
dom disorder were based on speckle patterns focused on 
the BEC • The large correlation length of the disorder 
created in this method could not cleanly reproduce the 
localization transition, although reduced diffusion was 
measured js], • The possible observation of this local- 
ization transition is further complicated by the fact that 
it is present only in three dimensional systems (in the 
thermodynamic limit). 

In this work we focus on a system consisting of a BEC 
in an optical lattice, to which a second optical lattice, 
with a wave vector incommensurate to the first, is added 
(13). This second lattice is weak compared to the first, 
and acts as a perturbation similar in certain aspects to 
random disorder. Specifically, for non-interacting par- 
ticles this model was considered in where it was 
shown that it exhibits a localization transition as a func- 
tion of the strength of the perturbation. Here we extend 
this model to the interacting case, and study the differ- 
ent phases as a function of both the interaction strength 
and the perturbation strength. A related setup, con- 
sisting of non-interacting atoms in two incommensurate 
lattices of equal strength, was studied in 11] . The prob- 
lem of bosons in two-color optical lattices was also dis- 
cussed in 12], in which commensurate lattices were stud- 
ied. The fundamental difference between commensurate 
and incommensurate lattices will be discussed later. Re- 



cently, the model studied here was realized experimen- 
tally ( 13]). We construct the phase diagram for incom- 
mensurate optical lattices, using an exact solution of the 
Bose-Hubbard Hamiltonian ([13) for small systems. We 
employ the superfiuid fraction and the determinant of 
the single-body density matrix [p) to distinguish between 
the phases [la]- It will be shown that for a mesoscopic 
system, the largest eigenvalue of p distinguishes between 
the phases as well. We then perform large scale mean 
field calculations and study the system in a realistic ex- 
perimental scenario. We show that the phase diagram 
exhibits the superfluid, Mott-insulator and Bose-glass 
phases, as well as effects of "disorder" induced order and 
delocalizing interactions. We also compare our system to 
the case of commensurate lattices, using both exact and 
mean field solutions of the Bose-Hubbard model, as well 
as GPE simulations. We then show that such incom- 
mensurability can be realized, and that our results are 
relevant for an actual experimental setup, taking into ac- 
count limitations of finite size and harmonic confinement. 

A system of bosons in a deep optical lattice in the tight 
binding regime is well described by the Bose-Hubbard 
Hamiltonian, given by |14| 
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where di^ and di are spatial creation and annihilation 
operators obeying bosonic commutation relations. J is 
the tunneling term between neighboring sites, Wi is the 
on-site energy at site i, and U is the on-site interaction 
energy. This model has been studied extensively [l^ . 
and is known to exhibit quantum phase transitions such 
as the superfluid - Mott insulator transition for Wi = 0. 

Following Aubry and Andre we first consider the 
homogeneous non-interacting limit, and introduce a sec- 
ond lattice with an incommensurate lattice vector as a 
perturbation. We choose the ratio of the wave vectors 
of the two lattices to be the golden ratio (3 — ^^^^.This 
perturbation alters the on-site energy term of Eq. to 
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Wi = Acos{2Tr(3i), giving 

iJ = - J ^ dj^di + A'^cos{2TTPi)di^d^. (2) 

<i,j> i 

Here A denotes the strength of the perturbation (or the 
depth of the incommensurate lattice). The addition of 
the second lattice alters the tunneling coefficient J as 
well, but we neglect this effect (this will be discussed be- 
low). This single-particle Hamiltonian was shown [lollr^ 
to have a localization transition in ID, with the critical 
point at A/ J = 2, by transforming to the momentum 
basis. It is important to note the similarity between this 
Hamiltonian and the Anderson Hamiltonian, with the in- 
commensurate lattice in the Aubry- Andre model playing 
the role of disorder in Anderson's model. However, an 
important difference between the models is that for dis- 
order, a delocalization-localization transition occurs only 
in three dimensions or more. In the Aubry- Andre Hamil- 
tonian such a transition can be seen already in one di- 
mension. 

We extend the Aubry- Andre model to include the ef- 
fect of interactions, by adding an on-site repulsion term 

H = —J aj^gj -|- A'y~^ cos{2tt [3i)di^ dj 



U 



(3) 



The Hamiltonian of eq. ^ can be solved by direct diago- 
nalization under periodic boundary conditions IT3] . How- 
ever, since the size of the matrix describing this many- 
body problem grows exponentially with the number of 
particles, we are limited to solving only small scale sys- 
tem, up to 9 particles in 9 sites. The properties of the 
system generally depend on the ratio of the number of 
particles to the number of sites, known as the filling fac- 
tor. In this work we focus on the case of unit filling factor, 
in which the Mott-insulator phase appears. 

In studying this system we use the superfluid fraction 
0, 0| and the determinant of the single-body density 
matrix [p) as the physical quantities which discriminate 
between the different phases 15]. The significance of p in 
describing the single-body properties of the system and 
the condensate fraction is detailed in ilSj. The quantity 
det{p), which is base-independent, serves as a measure 
of the amount of fragmentation in the system. If the 
system is completely fragmented, such that the occupa- 
tion of every single-body eigenfunction is non-zero, then 
det{p) ^ 0. Otherwise, the value of det{p) is zero. It can 
be seen that this quantity is related to the energy gap, 
which has been used as an order parameter in similar 
problems. This is a result of the fact that for complete 
fragmentation, the energy gap is non-zero (with a value 
of [/), and otherwise, since there are unoccupied states, is 
zero. However, det(p) is an advantageous criterion, since 
it emerges directly as a property of the ground-state. 

A widely used definition of the superfiuid fraction is 
through the response of the system to transverse drag- 



ging of its boundaries. The dragging is assumed to be 
an adiabatic acceleration process, in which the bound- 
aries reach a final, small velocity. This velocity can be 
described as a linear phase imprinted on the system, or 
a phase twist, and the superfiuidity as the phase rigidity 
of the system. It can therefore be calculated for a 
phase twist 9, with N particles and L sites through 



fs 



NJ 



02 



(4) 



In studying mesoscopic systems (as compared to taking 
the thermodynamic limit), det{p) can be replaced by the 
largest eigenvalue of p, which is the condensate fraction, 
in distinguishing between the phases. In the thermody- 
namic limit in the presence of interactions, the conden- 
sate fraction is zero both in the Bose-glass phase and 
in the Mott-insulator phase, and therefore cannot distin- 
guish between them. 
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FIG. 1: Det{p), superfluid fraction and condensate frac- 
tion as a function of interaction strength and perturbation 
strength, calculated by direct diagonalization of eq. ||3J for 
8 atoms in 8 sites. The first two quantities distinguish be- 
tween the superfluid phase {det{p) = 0, s/ 7^ 0), the Bose- 
glass phase {det{p) — 0,sf = 0) and the Mott-insulator phase 
{det{p) 7^ 0, s/ = 0). The resulting phase diagram exhibits 
disorder induced order for the regime of large interactions 
and delocahzing interactions in the regime of strong disorder 
(see text). Since the system is of finite size, the condensate 
fraction can also be used instead of det{p) to distinguish be- 
tween the phases, as it is zero in the Mott-insulator phase, 
and non-zero otherwise. 

In Fig. Qwe present the results of direct diagonaliza- 
tion of Hamiltonian ^ in one dimension. We show the 
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values of det{p), the superfluid fraction and the conden- 
sate fraction as functions of the interaction strength U 
and the strength of the perturbation (the ampUtude of 
the incommensurate lattice) A, both in units of the tun- 
neling energy J. Based on 15] we discriminate between 
the different phases of the system using det{p) and the 
superfluid fraction. Thus we note that for small inter- 
action strengths and small perturbations the system is 
in the superfluid phase, with det{p) = and the super- 
fluid fraction having a non-zero value. In this phase the 
mesoscopic condensate fraction (largest eigenvalue of p) 
is non-zero as well. By "mesoscopic" we mean a system 
of size TV, in which to every value of A/ J corresponds a fi- 
nite value of U / J, such that interactions weaker than this 
value give a condensate fraction cf >> 1/N. Keeping 
the strength of the perturbation small and increasing the 
strength of the interactions, the system passes into the 
Mott-insulator phase, where det{p) ^ and the super- 
fluid fraction is zero. Here the condensate fraction is also 
zero. For sufficiently strong perturbations, the system 
goes into the Bose-glass phase, in which the superfluid 
fraction is zero and det{p) is also zero. In this case the 
condensate fraction is positive [l5j |. Specifically, for zero 
interaction strength we recover the Aubry- Andre transi- 
tion from the superfiuid phase to the localized phase as 
a function of the perturbation strength, with the tran- 
sition point at A/ J — 2. Checking this result for other 
system sizes, wc find that the transition scales sharply 
to A/ J = 2, as opposed to Anderson's transition, which 
scales to A/ J = for ID 2, 15] . At the other extreme, 
for A = we recover the superfluid to Mott-insulator 
transition [l9|. It is important to note that det{p) and 
the superfluid fraction retain their distinction in the ther- 
modynamic limit, and are therefore valid order param- 
eters of the system. However, the condensate fraction 
scales to zero in the thermodynamic limit both for the 
Mott-insulator and for the interacting Bose-glass phase, 
and is therefore relevant only for mesoscopic systems. We 
have verified that indeed for larger systems c/ — > while 
sf —>■ const (for the normalization Tr{p) — N). 

The results of Fig. ^ show several similarities to 
the behavior of the Bose-Hubbard model with random 
quenched disorder. Most pronounced is the existence of 
the superfluid, Bose-glass and Mott-insulator phases and 
their relative position in phase space. Consequently, the 
effect of "disorder" induced order, i.e. the buildup of 
the superfluid fraction as a function of A in the strongly 
interacting regime, appears here. Also, the effect of delo- 
calizing interactions, which is the buildup of a superfluid 
response as a function of U in the regime of interme- 
diate perturbation strengths, is observed in this model. 
Finally, fragmentation of the condensate in the Bose- 
glass phase for a mesoscopic system, which occurs for the 
case of disord er lisll , is also found for an incommensurate 
perturbation 1201 ]. Fragmentation can be understood in 
the thermodynamic limit, by noting that above a certain 
strength of the perturbation, for every finite strength of 
interactions the lowest energy state will include a num- 



ber of fragments which scales with the size of the system 
[Tsl I2H . The transition from the Bose-glass phase to the 
Mott-insulator phase occurs when this fragmentation be- 
comes complete fragmentation, i.e. one particle per site 
at unit filling. 

Despite these similarities, the two forms of perturba- 
tions are inherently different, as indicated by the fact that 
in the Aubry- Andre model in ID the transition from su- 
perfluid to Bose-glass occurs at a finite value of the per- 
turbation, as opposed to the Anderson case, in which a 
transition exists only for three dimensions or more. An- 
other interesting feature, which is apparent in Fig. ^ 
are two stripes of a superfluid phase protruding from the 
central lobe of this phase. By studying the nature of the 
ground state of the system in this regime we find that 
these superfluid resonances are a consequence of ener- 
getic resonances between different many-body configura- 
tions. We have verified that such resonances exist for a 
variety of conditions, while their strength and position 
might vary. 

Thus far we have based our discussion on results ob- 
tained for a small system through direct diagonalization 
of Hamiltonian ijjji. To verify that our main results are 
not a consequence of finite size artifacts, we also per- 
form a large scale mean-field calculation based on the 
Gutzwiller ansatz '22] , adapted to the constraint of num- 
ber conservation. In this method we separate the many- 
body wavefunction to single-site wavefunctions 

iV')=nE/ni")' (5) 

i n 

where \n) are particle number states, and are the cor- 
responding amplitudes per site i. Minimizing the expec- 
tation value of the exact Hamiltonian for ]■;/') with 
respect to the /^'s, gives us the estimated ground state 
energy. In practice, in order to maintain a constant fill- 
ing factor, we need to add a chemical potential term 
—p'^i'Tii and to search for the minimizing /*'s in the 
restricted subspace where the total number of particles 
and the norms are conserved. We use the equivalence be- 
tween the Gutzwiller ansatz and the mean-field approach 
19, 23] to get around the non-linear minimization prob- 
lem. This is done by diagonalizing the mean field Hamil- 
tonian at each lattice point, and re-iterating this until a 
self consistent solution is reached, again while conserving 
both number and norm. 

In Fig. |2] wc show the results of the mean field cal- 
culation of the condensate and superfluid fractions as 
functions of both the interaction strength and the per- 
turbation strength for a larger system with 60 sites and 
particles. It should be noted that due to the decoupled 
nature of the Gutzwiller wavefunction, this method be- 
comes increasingly inaccurate as U/J lowers through the 
transition point and into the superfluid phase especially 
in ID. Moreover, it cannot reproduce the fragmentation 
of the condensate in the Bose-glass phase, which is seen in 
the exact diagonalization, and only sets an upper bound 
on the superfluid fraction (which is the average kinetic 
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FIG. 2: Superfluid fraction and condensate fraction as a func- 
tion of interaction strength and perturbation strength, cal- 
culated using the Gutzwiller based mean field approach for 
60 atoms and sites. The phase diagram can be qualitatively 
reconstructed, despite the limitations of mean field for low 
dimensional systems, near phase transitions and for the small 
interactions regime. The condensate fraction is seen to quali- 
tatively describe even the Bose glass phase for which it must 
approach zero at the thermodynamic limit, indicating the ex- 
istence of a mesoscopic regime. 
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FIG. 3: Superfluid fraction as a function of interaction 
strength and perturbation strength, calculated using exact di- 
agonalization (top), and the superfluid fraction as a function 
of A/ J for U/J = 1, calculated using the Gutzwiller based 
mean field approach (bottom), for commensurate (blue) and 
incommensurate (red) lattices. It can be seen that there is no 
dependence on A, and the superfluid fraction only exhibits the 
superfluid to Mott-insulator transition. The mean fleld cal- 
culation shows that the superfluidity remains unity as A/ J 
grows for commensurate lattices, while for incommensurate 
lattices it falls to zero. 



energy [22). For the same reason, resulting from the fact 
that this calculation does not conserve Tr{p), det{p) can- 
not be calculated here, and the condensate fraction must 
be used as an indication of the different phases. Never- 
theless, comparison between Figs. ^ and El reveals that 
the correct three phases of the system, as well as the su- 
perfluid phase resonances, are observed using mean field. 
The phase transitions themselves are smoothed out in 
the mean field approach, and occur at different values 
compared to the analytical result (for U = 0) and the 
exact solution (for U > 0). Although the condensate 
fraction must approach zero for a localized ground-state 
at the thermodynamic limit, we find that it is non-zero 
for systems with as many as 120 particles and sites (it 
does however become smaller for larger systems with the 
same parameters). 

The features of Hamiltonian ^ studied thus far are 
a direct consequence of the incommensurability of the 
two lattices. Similar properties were found in foi' 
commensurate lattices due to the finite size of the sys- 
tem examined therein. However, for thermodynamic and 
even mesoscopic systems, commensurate lattices exhibit- 
ing a periodic structure do not contain the phases found 
for incommensurate lattices. In Fig. 3a we plot the su- 
perfluid fraction as a function oiU / J and A/ J found us- 
ing an exact solution of Hamiltonian with (3 being a 
rational number, such that the two lattices are commen- 
surate. We also plot the superfluid fraction as a func- 
tion of A/ J for U/J = 1, found using our mean field 
approach, for both commensurate and incommensurate 



lattices. These calculations stress the difference between 
these two types of perturbations. Intuitively, incommen- 
surate lattices create a periodic structure, for which the 
unit cell is determined by the ratio between the wavevec- 
tors of the two lattices. Since the combined potential is 
still perfectly periodic, such a system should not exhibit a 
Bose-glass phase, effectively eliminating the effect of the 
perturbation strength A. In contrast, incommensurate 
lattices do not create a periodic structure, and therefore 
cannot be treated using Bloch theory. This fimdamen- 
tal difference gives rise to the localization phenomenon 
occurring as a function A, and the appearance of the 
Bose-glass phase. In accordance with this intuitive view, 
it is clear from Fig. 3a that the values of the superfluid 
fraction along the U/J axis are simply duplicated in the 
A/ J axis, indicating that the perfect periodicity of the 
system renders the A/ J axis redundant. This is in con- 
trast to the superfluid fraction given in Figs. 1 and 2, for 
which a dramatic change occurs in the A/ J direction, 
giving rise to localization and introducing the Bose-glass 
phase. 

We have also performed simulations of the Gross- 
Pitaevskii equation (GPE) in ID, in the regime of weak 
interactions (for which the GPE approximation is valid - 
i.e. for this mesoscopic system the interactions are weak 
enough such that fragmentation in the Bose-glass phase 
is negligible), and assuming "freezing out" of the radial 
dynamics. Using GPE, the external potentials of the in- 
commensurate lattices are used, and thus the values of J, 
U and A are derived and not introduced as parameters. 
The difference between commensurate and incommensu- 
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FIG. 4: Ground state real-space distribution for one lattice 
(blue solid), for a commensurate perturbation at A/ J ~ 50 
(red dashed-dotted), and for an incommensurate perturba- 
tion at A/ J ~ 50 (green dashed, normalized to 1/4 to fit 
scale), calculated using a ID GPE simulation. The localiza- 
tion is apparent for incommensurate lattices, as compared to 
commensurate lattices and to no perturbation. 



smaller than L. We have verified that this combined 
criterion can be fulfilled for realistic conditions, by cal- 
culating using both Gutzwiller mean-field and GPE the 
superfluid and condensate fractions aX U/J — 15 for a 
system with 60 sites, and lattices with wavelengths of 
78Anm and 1064nrn. The resulting behavior is the same 
as for pure incommensurability (using the golden ratio). 




rate lattices is manifest using GPE simulations as well. 
In Fig. 4 we plot the ground state spatial distribution for 
the case of one lattice (A/ J = 0), and for the case of two 
lattices (A/ J ~ 50) for both the commensurate and the 
incommensurate cases. It is clear that the ground-state 
of commensurate lattices is delocalized, with a structure 
resulting from the periodic unit cell of the combined lat- 
tices, while the ground-state of incommensurate lattices 
exhibits a well-localized single peak. 

Using our mean field method, we can further extend 
the Hamiltonian to better describe realistic experi- 
ments. The Hamiltonian Q can be approximately re- 
alized experimentally by loading a BEG into two opti- 
cal lattices with incommensurate wave vectors, one of 
which is weak compared to the other. This incommen- 
surate lattice may alter the tunneling energy J as well 
as the on-site energy U, an effect we have neglected so 
far. By using the Gaussian tight-binding approximation 
|14l[25l | we calculated the actual site-to-site tunneling co- 
efficients and the on-site interaction terms as a function 
of A. We have found that the variations in the interac- 
tion terms are completely negligible, and the tunneling 
energies vary by up to 0.2% for relatively weak perturba- 
tions (A < 30 J). In addition, GPE calculations for which 
the variations in J and U are taken into account inher- 
ently, show significant localization of the ground state as 
a function of the strength of the incommensurate lattice, 
which is apparent from both the spatial and momentum 
distributions (see Fig. 4). 

Experimentally, for finite size systems, the incommen- 
surability criterion states that the size of the unit cell 
of the combined lattices with momenta fci and /c2 must 
be larger than the size of the system. Moreover, the 
finite size L of the system dictates a momentum accu- 
racy Afc ^ fi/L which the system can distinguish. It 
is therefore necessary to add to the criterion a require- 
ment that for all the commensurate ratios p/q such that 
\p/q — ki/k2\ < Ak/k2, no p/q will have a unit cell 



FIG. 5: (a) Ground state spatial distribution calculated us- 
ing a 3D GPE simulation, for uj,, = 5kHz and /i — 2kHz, in 
the presence of one lattice (red) or incommensurate lattices 
(black). Compared to the truly ID geometry, the localiza- 
tion here is less pronounced, (b) Radial density profiles of 
the ground state for a high density regime (red) and a low 
density regime (blue). While the low density regime retains 
a gaussian profile, at high densities the profile assumes the 
Thomas-Fermi shape. 

An additional experimental consideration relates to the 
true dimensionality of the setup. We have thus far as- 
sumed a truly ID system, which requires tight confine- 
ment in the transverse dimensions such that they are 
"frozen out" . Naively, as long as the transversal trapping 
frequency ujr is larger than the chemical potential /z, the 
system can be treated as ID. However, since the ground- 
state spatial distribution is irregular in the presence of 
the incommensurate lattice, this argument is no longer 
true. In Fig. 5a we plot the ground state wavefunction 
calculated using a 3D GPE simulation (see, e.g., (23|)j as- 
suming radial confinement oi uj^ — ^ kHz and a chemical 
potential of /i = 2 kHz. Gompared to the ground state 
in the presence of incommenusrate lattices calculated as- 
suming a ID geometry shown in Fig. 4, the spatial dis- 
tribution here is clearly less localized. This suggests that 
despite the strong trapping in the transverse directions, 
the system does not retain its one-dimensional proper- 
ties in the presence of the perturbation. This results in 
a less localized state, with a smoothed transition com- 
pared to the ID case. To verify this we plot in Fig. 5b 
two radial profiles of the 3D ground-state, at longitudi- 
nal positions corresponding to high and low density cen- 
ters. It can be seen that while the low density regime re- 
tains a gaussian profile in agreement with the "transverse 
freezing" assumption, the high density regime displays 
a markedly Thomas-Fermi like profile, which indicates 
significant transverse dynamics. We therefore note that 
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standard magnetic or optical traps providing transverse 
trapping frequencies as high as a few kHz would not be 
able to create a truly ID setup, and would not exhibit 
a well pronounced localization peak. In order to achieve 
one-dimensional behavior in the presence of localization, 
an array of ID tubes created using a 2D optical lattice is 
needed. 




A/ J real space [lattice sites] 
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FIG. 6: (a) Superfluid fraction as a function of perturbation 
strength for constant interaction strength, calculated using 
the Gutzwiller based mean field approach for 80 sites and 
56 particles. U/J = 50, with a confining harmonic poten- 
tial, resulting in a smoothed and less pronounced transition 
from Mott-insulator to superfiuid. Nevertheless, the super- 
fluid fraction increases from approximately 0% to 10%, which 
may be detected experimentally. The dependence of the su- 
perfluid fraction on A/ J is sensitive to the specific form of 
the potential, but we have found that a response on the order 
of 10% is a robust result, (b) The particle number distri- 
bution in space for A/J = and A/ J = 9 reveals the so- 
called "wedding-cake" structure with Mott-insulator 
plateaus separated by superfluid regions. A small asymetric 
deformation of the structure can be observed as a result of 
the incommensurate perturbation. 

Finally, we include the inhomogeneity of a realistic sys- 
tem, by adding to the on-site energy the effect of an har- 
monic confining potential. The effect of the confining 
potential for a one-lattice system in the Mott-insulator 
regime (J7/J = 50) is to create a "wedding cake" struc- 
ture in space (in terms of the occupation number per 
site) [13, • By adding the incommensurate lattice as a 
perturbation, this structure is deformed, which results in 
an increase of the superfluid fraction of the system, ap- 



proximately from 0% to 10% (see Fig. 0). The superfluid 
fraction can be measured by dragging the lattices adia- 
batically, and through time of flight images distinguish- 
ing the atoms that remain stationary. The localization 
transition can be identified experimentally by observing 
a sharp increase in the three-body loss caused by the high 
density at localized regions 0. 

In conclusion, we have studied a BEC in an optical lat- 
tice, with a second, incommensurate optical lattice acting 
as a weak perturbation. We have calculated the phase di- 
agram of the system, using the superfluid fraction and the 
determinant of the single-body density matrix, and have 
recovered the three phases present in disordered mod- 
els - superfluid, Mott-insulator and Bose-glass. For a 
mesoscopic system, we have shown that the condensate 
fraction can also distinguish between the phases (instead 
of c?et(p)), and indicates fragmentation in the Bose-glass 
phase which increases with interaction strength, up to full 
fragmentation in the Mott-insulator phase. However, in 
the thermodynamic limit, the condensate fraction van- 
ishes in the Bose-glass phase for non-zero interactions, 
and therefore cannot act as an order parameter in that 
limit. Our solutions included both exact diagonalization 
of small systems and a new number conserving mean 
field method for large scale systems. We have comple- 
mented these calculations with ID and 3D GPE simu- 
lations. These allowed us to study realistic experimen- 
tal scenarios, including the effect of transverse trapping, 
varying tunneling energies and an harmonic confining po- 
tential. We have also used all three tools to demonstrate 
the fundamental difference between commensurate and 
incommensurate perturbations, as only the second ex- 
hibits localization. The similarities between this system 
in one dimension and disordered systems in three dimen- 
sions in terms of their phase diagrams, along with the fact 
that the perturbation we consider is created by introduc- 
ing a second optical lattice, suggest the possibility of an 
experimental realization of our proposed setup. Such a 
setup is simpler than a 3D one, and overcomes the corre- 
lation length limitations of disorder introduced through 
speckle patterns. 
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